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CONFLUENT VANDERMONDE MATRICES, DIVIDED 
DIFFERENCES, AND LAGRANGE-HERMITE INTERPOLATION 

OVER QUATERNIONS 

VLADIMIR BOLOTNIKOV 


Abstract. We introduce the notion of a confluent Vandermonde matrix with quater¬ 
nion entries and discuss its connection with Lagrange-Hermite interpolation over 
quaternions. Further results include the formula for the rank of a confluent Van¬ 
dermonde matrix, the representation formula for divided differences of quaternion 
polynomials and their extensions to the formal power series setting. 


1. Introduction 

The notion of the Vandermonde matrix arises naturally in the context of the La¬ 
grange interpolation problem when one seeks a complex polynomial taking prescribed 
values at given points. Confluent Vandermonde matrices come up once interpolation 
conditions on the derivatives of an unknown interpolant are also imposed; we refer to 
the survey m for confluent Vandermonde matrices and their applications. The study 
of Vandermonde matrices over division rings (with further specialization to the ring of 
quaternions) was initiated in [Tl]. In the follow-up paper [12] . the Vandermonde ma¬ 
trices were studied in the setting of a division ring K endowed with an endomorphism 
s and an s-derivation D, along with their interactions with skew polynomials from the 
Ore domain K[z,s,D\. The objective of this paper is to extend the results from [11] 
in a different direction: to introduce a meaningful notion of a confluent Vandermonde 
matrix (over quaternions only, for the sake of simplicity) and to discuss its connections 
with Lagrange-Hermite interpolation problem for quaternion polynomials. 

Let H denote the skew field of quaternions a = xo+ixi+ja^+kxs where xq, x±,X 2 , x% £ 
M and where i, j, k are the imaginary units commuting with M and satisfying i 2 = j 2 = 
k 2 = ijk = - 1 . For a £ HI as above, its real and imaginary parts, the quaternion 
conjugate and the absolute value are defined as Re(a) = xq , Im(a) = ix\ + jx 2 + ka^, 
a = Re(a)— Im(a) and \a\ 2 = oia = |Re(a)| 2 + |Im(a)| 2 , respectively. Two quaternions 
a and ft are called equivalent (conjugate to each other) if a = h~ l j3h for some nonzero 
h £ H; in notation, a ~ f3. It turns out that 

a ~ fd if and only if Re(a) = Re(/3) and |a| = \fl\, (1-1) 

so that the conjugacy class of a given a £ HI form a 2-sphere (of radius |Im(a)| around 
Re(cc)) which will be denoted by [cc]. It is clear that [cc] = {a} if and only if a £ M. 
For an a £ H\R, we denote by C Q the plane spanned by a and 1, which alternatively 
can be characterized as the set of all quaternions commuting with a. Observe that 
C a O [a] = {a, a}. 


i 
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Definition 1.1. A finite ordered collection ol = [a\,. .. ,ot k ) is a spherical chain (of 
length k ) if 

ot\ ~ «2 ~ • • • ~ Ofc and ctj+i fi&j for j = 1,..., k — 1. (1-2) 


The latter notion is essentially non-commutative: a spherical chain a consisting of 
commuting elements is necessarily a subset of C Ql fl [or] C {or, or} which together 
with inequality in (jl.2l) implies that all elements in a are the same: 

a = (a, a,... , a), a G EL (1-3) 

The spherical chain (11.31) may be understood as an element a taken with multiplicity k. 
Observe that many (if not all) results concerning multiple zeros of complex polynomials 
or multiple eigenvalues of complex matrices extend to the quaternion setting almost 
literally and can be interpreted in terms of commutative chains of the form (11.31) . It 
turns out, however, that they may have meaningful and nontrivial extensions to more 
general (non-commutative) spherical chains (11.21) . The latter observation is partly 
justified by the existing results on the zero structure of quaternion polynomials [4 
0 (9] and canonical forms for quaternion matrices 0 m 116] ; some extra evidence 
(the structure of the confluent Vandermonde matrices and interpolation conditions of 
Lagrange-Hermite type which do not appear in the commutative case) will be given 
below. 

1.1. The basic case. We now briefly review the basic (non-confluent) results in 
the form convenient for further “confluent” extensions. Let H[^] denote the ring of 
polynomials in one formal variable z which commutes with quaternion coefficients. The 
ring operations in Elfz] are defined as in the commutative case, but as multiplication in 1HI 
is not commutative, multiplication in IHIfz] is not commutative either. A straightforward 
computation verifies that for any a£i and / G H[z], 

f(z) = / e » + (z-a) ■ ( L a f)(z ) = f er (a ) + ( R a f)(z) ■ (z - a), (1.4) 

where f ei (a) and f er (a) are respectively, left and right evaluation of / at a given by 

m mm 

f et -{a) = ^a k f k and f er {a) = f k a k if /(z) = ^V/j, (1-5) 

k =0 k =0 j =0 

and where L a f and R a f are polynomials of degree m — 1 given by 

m— 1 / m—k—1 \ 77i—l / m—k—1 \ 

(L a f)(z) = a3 fk+j+l)z k , {Raf){z)=J2[ fk+j+lO.Az k . 

k =0 ' j =0 ' k =0 ^ j =0 ' 

(1.6) 

Interpreting Bljz] as a vector space over H, we observe that the mappings / i-A L a f 
and / i—^ R a f define respectively the right linear operator L a and the left linear oper¬ 
ator R a (called in analogy to the complex case, the left and the right backward shift, 
respectively) acting on H[z]. 

With evaluations (11.51) in hands, we formulate the left Lagrange interpolation prob¬ 
lem in H[z]: given distinct “points” or, ... ,a n G H and target values c±,... ,c n G H, 
find a polynomial f G H[z] such that 

f ei (ai) = Ci for i = l,...,n, 


(1.7) 



3 


or the right interpolation problem: find a polynomial f £ H[z] such that 

f ev (a.i) = Ci for i = l,...,n. (1.8) 

Making use of (11.511 (with m — 1 instead of m) we can write conditions (11.711 , (11.81) in 
the form 


/o + OLifx + ■■■a™ 1 fm—i =Ci (i = 1 ,..., n), 

fo + f\OLi + • • • fm- i«r _1 = °i (i = 1, • • •, n), 


respectively, with unknown /q, ..., f m - 1 , or in the matrix form, as 

ViF = C, F T Vff = C T , 

where F T = [f 0 ... f m -i ], C T = [c\ ... c n ] and where 



'1 

a\ 

a\ . 

_ m— 1“ 

. a 1 


1 


2 

_ m— 1 

V 1 = 

v m 

0-2 

a 2 • 



1 


a 2 

a m ~ l 


(1.9) 


are the left and right Vandermonde matrices associated with the given aq ,..., a n £ HI. 
Recall that the rank of a quaternion matrix is defined as the dimension of the left linear 
span of its rows or equivalently (by (Hi Theorem 7]), as the dimension of the right span 
of its columns. In general, the ranks of a matrix can be different from the rank of its 
transpose. For a set A, we will write jj(A) for its cardinality. The following result is 
due to T.-Y. Lam mi: 


Theorem 1.2. Let Si,... ,Sf be all distinct conjugacy classes having non-empty inter¬ 
section with the set A = {aq ,... ,a n }, and let 

f i, if MS} n A) = l, .. 

K = ^ + ... + fH, where H = [ 2 , if fl(S-n A) > 2. (L10) 

Then rankV)^ = rankfY) = min(m, k) . In particular, the square matrix is invertible 
if and only if all elements in A are distinct and none three of them belong to the same 
conjugacy class. 

The result was established in mi in a more general setting of division rings with a 
fixed endomorphism and the integer k was identified with the minimally possible degree 
of a nonzero polynomial having left zeros at ai,...,a n . It follows from Theorem 11,21 
that if none three oi a\,... ,a n belong to the same conjugacy class, then the problems 
G2D and (|1.8I) have (unique) solutions of degree less than n for any choice of off s. 
Otherwise, the problems may have no solutions, which is indicated by the rnext result 
( 0 ). 

Theorem 1.3. For f £ lH[z] and three distinct equivalent elements a 1 , 02,03 € H, 

f e£ (03) =(a 3 - a 2 )(ai - a 2 )~ 1 f ei (a 1 ) + (oq - a 3 )(ai - 02)^f e£ (a 2 ), (1.11) 

f er (a 3 ) =f er (a 1 )(a 1 - a 2 )~ l (a 3 - a 2 ) + f er (a 2 ){ai - a 2 ) _1 («i - a 3 ). (1.12) 
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Indeed, in order the problem (11.711 or (jl. 8 j) to have a solution, the target values 
assigned to any triple of equivalent nodes must satisfy certain consistency conditions 
pointed out in Theorem 11.31 If the problem is consistent, then it admits a (unique) 
solution of degree less than rank 17^; see [3j for details. 


1.2. Main results. Left divided differences of a given polynomial based on the 
ordered collection a = (or, ..., a &) of nodes are defined in terms of left backward shift 
operators (USD by [ay,f]e = / e *(oi) and 


[ai,...,a j ;f] £ = {L olj _ 1 ---L Ql f) ei (aj) for j> 2. 

The right divided differences are defined by similar formulas (with right evaluations 
instead of the left and with R a instead of L a \ see (HID). As we will see in Section 
2 , quaternionic divided differences are quite different from their complex prototypes, 
unless the nodes ctj’s commute, i.e., ct? £ C ai for i = 1,..., k. Divided differences based 
on a spherical chain are of special interest. Let us introduce notation 


A i(or,f) 


[«i;/k 

[ai,a 2 ;/b 

[ai,a 2 ,■■.,«&;/]£ 


a = (ai,... ,a k ), f eU[z] 


(1.13) 


for the column of left divided differences of a given polynomial / based on the spherical 
chain a. The next theorem shows that the columns A i{ocy /) and Ai(a 2 ', f ) associated 
with the chains of the same length k from the same conjugacy class Sci and with 
distinct leftmost entries determine all left divided differences of / of order up to k based 
on any chain 0:3 C S. 


Theorem 1.4. Let ex? = (a* 1 ,..., a, &) (i = 1,2,3) be three spherical chains of the 
same length and in the same conjugacy class Sci and let 04,1 f 0 : 2 , 1 - Then there 
exist k x k matrices A and B such that 

^e{^ 3 ‘,f) = AAi(oLyf) + BAi(oL 2 -,f) for any f £ M[z}. (1.14) 


The proof (along with explicit formulas for A and B) will be given in Section 5. 


Definition 1.5. Given a spherical chain a = ( 04 ,..., o^), we define the left confluent 
Vandermonde matrix V^(cx) by 

(1.15) 

(1.16) 

(1.17) 


^m(^) [jo^i ? 0(2) • • • j ^ ^ 

The left confluent Vandermonde matrix based on n spherical chains 

LX? — (o?, 1 , • • • , ), i — 1 , ... , Tfc, 

is given by 


V£(ai,...,a n ) = 


V*{<* 1 ) 




where the matrices V^iaf) are defined via formula (11.151) . Similarly, the right confluent 
Vandermonde matrix based on spherical chains (jl.161) is given by 

V£(a 1 ,...,a n ) = [V^(ai) ... V^(a n )) , (1.18) 
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where the m x k block corresponding to a sole chain a = (ai,..., a^) is defined as 

Ki( a ) = [[^; «1, «2, ■ ■ ■ , • (i- 19 ) 


Vaguely speaking, the confluent Vandermonde matrix should be defined so that it 
will be non-singular in situations where the usual Vandermonde matrix is singular, that 
is, according to Theorem 11.21 if 


(1) a n = a.i 2 and/or (2) a. n ~ a.i 2 ~ aj 3 . 


It turns out that the matrices (11.171) . (I1.18P do the job. As we will see in Section 
4, these matrices arise in the context of Lagrange-Hermite interpolation, like regular 
Vandermonde matrices do within Lagrange interpolation problems (11.81) . (11.91) . In case 
ki = 1 for i = 1,...,n, the matrices (11.171) . (11.181) amount to V^ n and V/j in (|1.9D . If 
the spherical chain a is of the form (11.31) . then V^(a) and V^a) take the form 



”1 a a 2 . a m 1 

0 1 2a . (m—l)a m-2 



1 ka ... 


(m— 1 )! ^m — k 
(n-l)! " 



( 1 . 20 ) 


In the contrast to the basic case (HD, left and right confluent Vandermonde matrices 
based on the same noncommutative spherical chains in general are not transposes of 
each other (see Example 12.61 belowh However, they are related as is indicated in the 
following result (see Section 2.2 below for the proof). 


Lemma 1.6. Let V)^(q;i, ..., a n ) and V^(a i,... ,a n ) be the left and the right confluent 
Vandermonde matrices based on on spherical chains ai = (a^i,..., and ai := 

(a^i,..., respectively. Then 

^(ar, ■ ■ ■, a n ) = (^(a b ..., a n ))*. (1.21) 


To present the confluent version of Theorem 11.21 we need the analogs of the integers 
pfl s in (11.101) . If the conjugacy class S C HI contains only one spherical chain ccj from 
(11.161) . we let p(S) = ki. If S contains d > 2 spherical chains (|1.16l) . we pick the longest 
chain ct, = (a^i..., C S, and for any other chain aj = (a^i..., OLj^d) C S, 

define the integer 


0, if otj t i ^4 Q^i, 

max{r : aj£ = oti£ (1 < ^ < r)}, if aj t i = cc^i. 


( 1 . 22 ) 


We then let 


/r(5) = ki + maxjfcj — Uj}. (1.23) 

Note that if there are several chains of the maximal length, then the values of Vj in 
(11.221) depend on which one of the longest chains has been chosen for the comparison. 
However, the integer (11.231) is independent of this choice (see Proposition 13.51 below). 
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Theorem 1.7. Let V^(qi, ..., a n ) and V^(qi, ..., «„) 6e t/ie confluent Vandermonde 
matrices based on spherical chains (11.161) . To each conjugacy class Sj containing at least 
one of these chains, assign the integer p(Sj) as in (ll.23j) . fjl.22[) . and let k = )TL p(Sj)- 
Then 

rank i,..., a n ) = min(m, k) = rank V^(qi, ..., a n ). (1-24) 


Remark 1.8. The second equality in (|1.24l) is a consequence of the first, by Lemma h. 6 1 
Indeed, upon applying formulas (jl ,22[) to the chains a i,..., a n we come up with the 
same integer k as for the original chains or, ..., a. n . Assuming that the first equality in 
(11.241) holds true, we conclude that rank V^(qi, ..., a n ) = rankI^(cEi,..., a n ). Since 
rank A = rank A* for any matrix A over H, the second equality in (11.241) follows by 

(fOIT) . 

The outline of the paper is as follows. In Section 2, we present confluent Van¬ 
dermonde matrices as unique solutions of certain Stein equations and prove Lemma 
11.61 In Section 3, we recall some basic facts on indecomposable polynomials and show 
that the integer k in Theorem 11.71 can be alternatively introduced as the degree of the 
least right common multiple of polynomials P ai ,... , P an associated with the spheri¬ 
cal chains oq, ..., a n by formula (13.41) . In Section 4, we demonstrate how confluent 
Vandermonde matrices arise in the context of a Lagrange-Hermite type interpolation 
problem and present necessary and sufficient condition for a square confluent Van¬ 
dermonde matrix to be invertible. The proofs of Theorems 11.41 and 11.71 are given in 
Section 5. Several extensions of these results to the setting of formal power series over 
quaternions are presented in the concluding Section 6. 


2. Preliminaries 


In the complex setting, divided differences are defined in terms of the operator L a : 
f{z ) —> acting as a backward shift on the sequence of Taylor coefficients of / 

at a given point a £ C. Evaluation formulas for L a f are: 


Laf(fl) = 


/(£)-/(«) 

(3—a 

/'(«) 


if fl / a, 
if fl = a. 


( 2 . 1 ) 


Furthermore, the operators L a and Lp commute and satisfy the Hilbert identity Lp — 
La = (ft — a)LpL a due to which L ari _ 1 ... L a2 L ai f(a n ) is the n-th divided difference of 
/ at the nodes or, ..., a n which in turn, is an ingredient of the Newton’s interpolation 
formula. 


In the noncommutative quaternionic setting we distinguish the left and the right 
backward shifts (11.61) . unless a is real in which case L a = R a . If we denote by /^ the 
k -th formal derivative of / G H[z], then a straightforward verification shows that for 
any fixed a£i, 


deg(/) 


/- E 


k =0 


(/ (fc) )^(a) 

k\ 


deg(/) 


v- (f (k) ) er (u)„k 

/ . f.i Pai 


k =0 


k\ 


p a (z ) := z — a. 


( 2 . 2 ) 
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In terms of the latter Taylor expansions, the operators (| 1. 6 11 take the form 

deg(/)-l /j-Cfc+mp.deg(/)-l (fc+1) r ( a ) 

R-af = ^2 


L a f= Y. 


pi 


k=0 


(/( fc+ 1 )) e *(a) 
(* + !)! 


k =0 


(k + 1)! 


k 

on 


which justifies the “backward shift” terminology. The operators L a and Lp are essen¬ 
tially non-commuting. 


Proposition 2.1. L a Lp = LpL a if and only if a(3 = (3a, in which case 

La — Lp = (a — (3)L a Lp. (2-3) 

Indeed, if a(3 = (3a, then the asserted equalities are verified as in the complex case. On 
the other hand, since (L a Lp — LpL a )(z A ) = (3a — a(3 (by the first formula in (11.61) 1. 
commutation equality L a Lg — LpL a implies (3a = a(3. A similar statement holds for 
right backward shifts. 


2.1. Divided differences. Given a polynomial / € H[^], the successive application of 
formula (11.41) to elements a\,... ,a n £i and polynomials /, L ai f, L a2 L ai f,... leads 
us to the representation 

71—1 

/ = r*{a r) + £ Pai ... Pak ■ ( L ak ■ • • L ai fr*(a k+1 ) 

k= 1 

+ Pai ■ ■ ■ Pa n ' (La n ■ ■ ■ L ai f), (2.4) 

which, being the (left) quaternionic analog of the Newton interpolation formula, sug¬ 
gests to define quaternionic left divided differences 

[ai;/k = f e£ {ai), [ai,...,a k -,f] £ = (L ctk _ 1 ---L ai f) e£ (a k ) for k > 1. (2.5) 

Similarly, the formula based on the successive application of the second representation 
in (fT~41) . 

71—1 

/ = f £r {ai) + ^2(Ra k ■ ■ ■ R a J) er (a k + 1) ■ Pa k - Pai 

k=1 

+ (Ran ' ‘ ‘ Raif) ' Pa n • • ■ Pai (2-6) 

suggests to introduce quaternionic right divided differences by 


[/;aq]r = f er (&i), [f-,ai,...,a k } r = (R ak _ 1 ---R a J) er (a k ) for k > 1. (2.7) 


Letting aj = a for j = 1,..., n in (12.41) . (12.61) and comparing the obtained representa¬ 
tions with ( 12 . 21 ) we conclude that 


[ a,. .j , cr , 
(/c+1) times 


(f {k) ) ei (g) 

k\ 



(fc+1) times 


(/ (fc) ) e ”(«) 

k\ 


( 2 . 8 ) 


for all k > 0. The latter formulas justify equalities (11.131) . 

The difference between the complex and quaternionic settings becomes transparent 
even in the case where k = 2. It is not hard to show that if a 2 aq, then 

[oq, a 2 ; f]e = (a 2 - aq) _1 (/ e *(a 2 ) - f e *(a i)), 


(2.9) 
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where a 2 = (a 2 ~ Sq) - 1 02(02 — ar). The formula (12.91) is similar to its complex 
counterpart ( 12 . 11 ) . but the element 02 is replaced by the equivalent element a 2 which 
is equal to 02 if and only if oq and «2 commute. If aq ~ 02 / aq , then (as we will see 
in Example 15.71 below) the formula for [oq,a;2; f]e. invokes not only the values of / at 
aq and 02 , but also the value of f: 

[a 1 ,a 2 ;f]e = (o 2 - a 2 ) _1 (/ e£ (a 2 ) - / e£ (oq) + («2 - ai)/ /e£ (aq))). (2.10) 

Thus, the divided differences based on a spherical chain (different from that in (12.81) 1 
is the object which does not appear in the commutative setting. 

Lemma 2.2. For any h £ H[z], x £ M and aq, 02 ,..., ajt £ H, 

[aq, • • •, otk 1 — (a.k x ) [ai,..., a^, L X K\^ T [ai,..., otk— 1 , . (2.11) 


Proof: Proposition 12.11 applies to a = ak-i and /3 = x (since x is real), and the 
formula (12.31) implies 

La k _! = (I + {oik -1 — x)L ak _ 1 )L x . 

Applying formula (11.41) to / = L ak _ 2 ■ ■ ■ L ai L x h and a = ak -1 leads us to 
^~‘ a k —2 "" ’ ^ j a\L x h {z oik—\)L ak _ 1 L ak _ 2 • ■ ■ L ai L x h 
= {La k _ 2 ‘ ‘ ‘ L ai L x h) 1 {oik— 1 ) = [oq j ■ ■ ■, oik— lj L X K\^ . 

Making use of the two last relations and taking into account that L x L a:j = L a .L x (by 
Proposition ^. II) . we have 

L ak _ 1 L otk _ 2 ■ ■ ■ L ai h {z x)L ak _ 1 L olk _ 2 • • • L ai L x h 

— (-1 T {oik—i x)L ak _^L x L ak _ 2 • • • L ai h 

{z oik—i)L ak _ 1 L olk _ 2 • • • L ai L x h {oik—\ x)L ak _ 1 L olk _ 2 • • • L ai L x h 
— L x L ak _ 2 • • ■ L ai h L ak _ 2 • • • L ai L x h T [aq,..., oik—i , E x h\g 
= [ai,..., oik —1 ; L x h\g • 

Evaluating the latter equality at z = ak on the left gives, in view of ([2.51) . 

[aq: • • • 1 Ok} h] g {oik *^) [ai, • • •, oik 1 L x [ai,..., a^— \ , L x 

which is equivalent to (12.111) . □ 


2.2. Explicit formulas. In this sction we present the explicit formulas for matrices 
V^{a) and V^{a) defined in (11.151) . ([1.191) . Here we will not assume that the collection 
of nodes a = (aq,..., a*,) is a spherical chain. Given a = (ai,..., a*,), we let 


0 

0 

S 



1 02 0 


1 



0 

0 1 : 

and Ek = 


: 0 


0 

0 ... 0 1 afc 


L - 1 


and we let Vij denote the (i. j)-entry of the matrix V^{cx)\ 

Vm { a ) = > v i,j : = [<Xl,<*2, Z 3 ] t . 


( 2 . 12 ) 
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Since deg (L a f) = deg(/) — 1 for each polynomial / of positive degree and since the 
backward shift of a rnonic polynomial is again monic (or identical zero), it follows 
that V^(ck) is upper triangular with all diagonal entries equal one. In particular, the 
leftmost column V\ of V^a) equals E k . Applying equality (|2.11|) to x = 0 and h = z 3 
and taking into account that Lqz 3 = z 3-1 , we get 

[«i,.. .,a k \z J ] l = a k [«i,.. ,,a k ;z 3 ~ 1 ] £ + [«i,..., a fe _i; z J ~ l ] £ , (2.13) 

which implies the recursion 

Vij = diVij-i + (■ i,j > 2) (2.14) 

for the entries of V^(a). Observe that 

Oj = K;^] = «i = ai^ij-i, 

which together with (I2.14|) imply that the consecutive columns in V^ot) are related by 
Vj = JcxYj-x (j = 2,... ,m). Since V\ = E k , the latter recursion gives V) = 1 £'fc. 

Similarly, using relation 

[^;ai,... ,a k \ r = a k [z 3 ~ l \ «i,..., a k ] r a k + [V -1 ; ai,..., a fc _i] r , 

the “right” counter-part of (12.131) . one can show that the j-th row in the matrix (11.191) 
equals Ej(JI .) 3 ^ 1 . We thus arrive at the following result. 


Lemma 2.3. Given a = (ai,..., a k ), let V~flpL) and Vffla) be defined as in (11.151) . 
(|1.19l) (confluent Vandermonde matrices if a. is a spherical chain). Then 


Vi{cx) = [E k J a E k ... JZ~ l E k \ , t£(a) = 

where J a and E k are given in (12.121) . 


El 

W t~X 


E k 3 ol 

Ei(j^r~\ 


, (2.15) 


Remark 2.4. If we let F m be the m x m lower triangular Jordan block with zero 
entries on the main diagonal (F m = [<5j_ij]™ =1 , where dij is the Kronecker symbol), 
then we can rewrite formulas (12.15|) as 

771—1 771—1 

V*{a) = E J a E k El{Fly, V^(a) = E F^eJ (J^) 3 , (2.16) 

3 =0 3 =0 


so that V^(a) and Vfifioc) are unique solutions to the respective 






V*(a) - J a V*(a)FZ = E k El, K)(a) - F m V^{a)J^ = E m Ej. (2.17) 


Proof: Indeed, if A is a solution to the Stein equation X = J^XF^ + E k Ej n , we 
can iterate this equation to represent X in the form 

OO 

X = E J a E k El(FZ) j 

3=0 
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and then to observe that since = 0, the formula for X is the same as that for 

V^(a) in (12.161) . The fact that Vffa) is a unique solution to the second equation in 
(12.171) is justified similarly. □ 


Remark 2.5. Since left evaluation functionals and left backward shift operators are 
right linear, the formula (11.131) defines a right linear operator / i— > A^ (a; /) acting from 
IHIfz] into M k . Again, we do not assume that a = a±,... ,a k ) is a spherical chain. By 
linearity, we have 


m— 1 

A*(a;/) = J 2 

3=0 


[ar,z j ]i 

[ai,a2;z J ’]f 


m—1 

fj, if f(z) = ^2 

3=0 


_ [a\,a2,—,a k -,z J ]i . 

which on account of definition (|1.15l) and formula (12.151) . can be written as 


m—1 


A t(cr,f) = ' 52 jiE k f j = V*(ot) 

3=0 



(2.18) 


Proof of Lemma II. 61 Due to definitions (11.171) . (11.181) . it suffices to verify (11.21[) 
for the left and right confluent Vandermonde matrices based on a sole chain, i.e., to 
verify that 

V£(a) = 0C(a))* for a = ( or , . .., a k ). (2.19) 

The latter matrices are unique solutions to equations (12.171) . Replacing a by a and 
taking adjoints in the second equation in (|2.17|) we get (since (l7j)* = J a ) 

CC(a))* - J a (VZ(a))*FZ = E k El. 

Therefore, V^(a) and (V^(a))* satisfy the same first equation in (12.171) . which has a 
unique solution so that V^(a) = (V^j(«))*. □ 

Example 2.6. Let us apply formulas (12.151) to the spherical chain a = ( a,/3 ) (a ~ 
/3 a) to get 



T 

a 

a 2 


a 3 


a m ~ l 


vi(cc) = 

0 

1 

(3 - h Q; 

1 3 ’ 2 

+ / 3 ct + a 2 

m— 2 

■ E 

p3 a ™-i~ 2 

, ( 2 . 20 ) 







3=0 




T 

a 

a 2 


a 3 


a m ~ 1 


(C(a)) T = 

0 

1 

(3 Qi. 

fl 2 

+ 0.(3 + Q 2 

m— 2 

■ E 

a rn-j-2pj 








3=0 




Examining the (2,4)-entries shows that V/(a) V 4 r (ct) T unless a/3 = /3a which is 

possible only if (3 = a. On the other hand, it is readily seen that V 4 (a) = 

(which agrees with Lemma 11.61) . 

2.3. Similarities and distinctions with the commutative setting. In the com¬ 
mutative setting, confluent Vandermonde matrices can be obtained from the basic ones 
by certain limit procedure, sometimes called the confluence of one row into another. 
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We now examine this procedure on 2 xm quaternionic left Vandermonde matrices. Ac¬ 
cording to the commutative recipe, we start with the left Vandermonde matrix based 
on the elements a and a + e and then pass to the limit as e —> 0 in the following 
product: 


lim 
£—>0 


( 2 . 21 ) 


1 a a 2 ... a m 1 

1 a + e {a + e) 2 ... (a + e) m_1 
The two leftmost elements in the bottom row are 0 and 1. To compute other entries, 
we observe that 


1 

0 ' 


-e- 1 




k -1 

£ —1 ^(a + e) k — a k ^j = £ _1 a^ ea k ~^~ l + e _1 i?(a, e) (2.22) 

3 =0 


for each k > 2, where R(a, e) is the sum of all words of length k in letters a, e and 
containing at least two letters e. Then it follows by the triangle inequality that 



as e —1 0. 


We now conclude from (12.221) that the limit (12.211) exists (for m > 3) if and only if the 

limit 3 = lim £ _1 oi£ exists, that is, never , unless a £ M (in which case the confluent 
£—^0 

matrix takes the form (11.201) 1. However, we may adjust (I2.2ip by taking restricted 
limits (i.e., by sending e to the origin along continuous curves) calling any matrix 
arising in this way a confluent Vandermonde matrix. If such a restricted limit 3 exists, 
it necessarily belongs to the conjugacy class [a] and the corresponding restricted limit 
(|2.21l) amounts to the matrix V^a.) from (|2.20l) . To show that any matrix of the form 
p.201) (that is, with any 3 £ [a]) arises in this way, it suffices to represent 3 i n the 
form 3 = hah,- 1 for some h ^ 0 and evaluate the limit (J2.21I) as e —)■ 0 along the line 
= {fh : r £ M}. 


3. Indecomposable polynomials 

In this section, we record several basic facts on quaternion polynomials needed in 
the sequel. As the division algorithm holds in H[z] on either side, any (left or right) 
ideal in H[z] is principal. We will use notation ( h) r := {hq : q £ H[z]} for the right 
ideal generated by h, and we will write / = g (mod r h) in case (/ — g) £ ( h) r . Similar 
notations will be used for left ideals. Maximal ideals in H[z] are generated by linear 
polynomials p a (z) := z — a (a £ H): it follows from (11.41) that 

/ £ <P Q )r & f ei {a) = 0 and / £ {pp)i 44 / e -(/3) = 0. (3.1) 

In the latter cases we say that a and 3 are respectively, left and right zeros of /. As 

was shown in M, any (monic) polynomial p £ H[z] of degree deg(p) = k > 1 has a left 

(a right) zero, which along with ()3.1j) implies that / can be factored into the product 
of linear factors 

p = p ai ...p ak , p a .(z) = z-aj. (3.2) 

Equivalences (13.11) are particular cases of more general statements 
/ e {pen ■ • • PaJ r 44 [«i,.. .,aj\ f] £ = 0, 

/ € (p ai • • • p a Je 44 [/; ai,..., ay] r = 0, 

which in turn, follow from (12.4(1 — (12.7(1 . 


for j = 1,..., k, 


(3.3) 
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Given polynomials f,g £ Hfz], their least right common multiple h = lrcm (f,g) 
and their least left common multiple h = llcm(/, g) are defined as (unique) monic 
polynomials such that 

Wr = (/) r n(s) r , (h)t = {f)t n (g)t. 

Following [ T5] . let us say that a polynomial / is indecomposable if it cannot be rep¬ 
resented as the lrcm of its proper left (equivalently, as the 11cm of its proper right) 
divisors. The latter means that the ideal (h) r is irreducible in the sense that it is not 
contained into two distinct proper right ideals in Hfz]. Various characterizations of 
indecomposable polynomials are listed below (see e.g., [4] for the proof). 

Theorem 3.1. Let p £ H[;s] be factored as in (13.21) . The following are equivalent: 

(1) The ideal (p) r is irreducible, i.e., p is indecomposable. 

(2) The ideal (p)i is irreducible. 

(3) a = (ai,..., otk) is a spherical chain. 

(4) ot\ is the only left zero of p. 

(5) is the only right zero of p. 

(6) (13.21) is a unique factorization of p into the product of linear factors. 

The latter theorem establishes a one-to-one correspondence 

a = (ai,... ,a k ) e-> P a = p ai .. ,p ak (3.4) 

between spherical chains and monic indecomposable polynomials. 

Given f,g£ H[;z], their greatest left common divisor glcd(/, g) is defined as a monic 
polynomial d of the highest possible degree such that f = df and g = dg for some 
f,g£ H[;z] or equivalently, as a monic generator of the right ideal (/) r + (g) r . 

Remark 3.2. It follows from the equivalence (3) <(=> (6) in Theorem 13.II that given two 
spherical chains on = (a«,i,..., a, ^) (i = 1, 2) belonging to the same conjugacy class 
and sharing v leftmost elements, the glcd of the associated polynomials P ai and P a2 
(see (1231)) equals p Q] • • • p a „. 

Let us recall that the characteristic polynomial of a nontrivial conjugacy class S C H 
is defined by 

X$(z) = (z — a)(z — a) = z 2 — z(a + a) + |a| 2 , (3-5) 

where a is any element in S; it follows from characterization (11.11) that formula (|3.5j) 
does not depend on the choice of a € S'. Since Xs is the polynomial of the minimally 
possible degree such that its zero set (left and right, as X £ R[z]) coincides with S, it is 
also called the minimal polynomial of S. We now recall a result from [4] concerning least 
common multiples of indecomposable polynomials having zeros in the same conjugacy 
class. 

Lemma 3.3. Given indecomposable polynomials 

Pa.i Pcti t i • • • Pcti'b.. i (ch,l) ■ ■ ■ ) ) i (3-6) 

based on spherical chains cx\,... ay in the same conjugacy class S such that deg(P Q , 1 ) > 
deg(Pa, 2 ) > ... > deg(P ad ) (i.e., h > k 2 > ... > k d ), let 

Polj =Pjhj, where pj = glcd(P a!j ., P ai ) for j = 2,...,d, 


(3.7) 
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and let m = max deg(hj). Then 
2<j<d 


lrcm (P ai ,...,P ad ) 


ym 
*S > 

ym 

Pcn,iPai,2 ' 


P, 


C^k-i —r, 


if m = ki, 
if m < k\. 


(3.8) 


Corollary 3.4. If on = (cc^i,..., (i = 1,2) are two spherical chains in the same 
conjugacy class and i/ai,i ^ 02 , 1 , then deg(lrcm(P ai , P a2 )) = deg(P ai ) + deg(P a2 ). 

Since all polynomials (13.61) are indecomposable, it follows by Remark [3 .2 1 that deg (pj) = 
Uj, where Vj is the integer defined in (11.221) (with i = 1 since in Lemma 13.31 we assumed 
that aq is the longest chain). Therefore, the integer m in (13.81) equals max (kj — Vj ). 

By formula (13.81) , 

deg(lrcm(P ai ,... , P ad )) = k\ + m. 

Comparing the latter equality with (|1.23l) leads us to the following 

Proposition 3.5. The integer p,(S) defined in (11.231) . (11.221) is equal to deg(lrcm(P ai : 
on C 5)). 

We next observe that if f\,...,f n € BI[ 2 ] are “spherically coprime” (i.e., no ze¬ 
ros of fi and fj belong to the same conjugacy class), then deg(lrcm(/i,..., f n )) = 
deg (fj). Combining this observation with respectively, Corollary 13.41 and Propo¬ 
sition [33] leads us to the following conclusions. 

Remark 3.6. If the leftmost elements agi,..., a n ^\ in the chains (|1.16|) are all distinct 
and none three of them belong to the same conjugacy class, then 

n n 

deg(lrcm(P Ql ,..., P a J) = ^ deg (P aj ) = ^ k r 

3 =1 3 =1 


Remark 3.7. The integer k in Theorem 11.71 equals deg(lrcm(P Ql , ..., P an ))- 
4. Lagrange-Hermite interpolation 

In Section 1, we recalled how Vandermonde matrices arise in the context of the 
Lagrange interpolation problem. In this section we formulate and solve a more gen¬ 
eral Lagrange-Hermite problem and will see (in formulas (|4.4|) and (14.61) below) how 
confluent Vandermonde matrices arise in this more general setting. 

Due to equivalence (13.11) . the typical interpolation condition f ee (a ) = c in the left 
Lagrange problem gzd can be written in the form / = c(mod r p Q ,) where the target 
value c is understood as an element of H[z] of degree zero. This form suggests to 
consider a more general condition 

/ = h (mod r p) (4.1) 

where p is a given indecomposable monic polynomial and h is a given polynomial of 
degree deg (h) < deg (p) = k. By Proposition 12.11 we may take p in the form (|3.2|) 
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where a. = (or, ..., otk) is a spherical chain, and we may apply representation ([2.411 and 
notation ra to write 

k -1 

h = hi + ^ p ai ■ ■ ■ p a . • fej+i, where hj = [or,..., ay; /i]*- 
i=i 

Since any / £ H[z] can be represented as 
k -1 

/ = K; /k + • • • Paj ' [° lj • • • > "d •fk + pq for some Q e 

J=1 

we conclude that condition (14.111 can be equivalently written as 

[a 1 ,...,a j ;f]i = h j for j = l,...,k, (4.2) 

or in the vector form, upon making use of notation (11.131) . as 


A e{a;f) = Cob <j<khj. (4.3) 

Upon invoking formula ([2.1811 for A^ct;/), we write interpolation condition (14.31) in 
terms of coefficients of the unknown / as 




fo 


" hi “ 

_ fm—1 _ 


_ hk _ 


(4.4) 


Similarly, given an indecomposable polynomial p = p Qfc ... p ai (observe that according 
to definition if a = (ai,...,afe) is a spherical chain, then the reversed tuple 

(afc,...,ai) is also a spherical chain) and given h £ H[z] (deg(/r) < k), the “right” 
analog of condition (14.11) is equivalent to the matrix equation 


f = h (mod e p) [f 0 ... fm-i] V£(a) 



where hj = [hr, or,..., aj\ r for j = 1,..., k. 


Remark 4.1. For the spherical chain (|1.3I) . interpolation conditions (14.2p prescribe the 
values of (a) for j = 0,..., k — 1 (see formulas ([2.811 1 bringing up therefore, the 

literal quaternionic analog of the well-known Lagrange-Hermite interpolation problem. 
In the general case, as is indicated in (12. 10I1 . conditions (14.21) prescribe certain combi¬ 
nations of the values of / and its derivatives at the elements of the chain a. Since the 
ideal (p) r is irreducible, the values of / and its derivatives cannot be separated, so we 
indeed have an interpolation problem which does not appear in the commutative case. 


Left Lagrange-Hermite interpolation problem: given n spherical chains (11.161) 
and given elements Cij € H, find a polynomial f £ H[,z] such that 


i ki■ 


■ ■ ■ ) a i,j j f]i — c i,j f or * — \, ... ,1V, j — 1, 

Similarly to (14.41) . we can write conditions ([4.511 in the matrix form as 


f / m(^-l? • • • ? (X-n) 



' fo ' 


'Ci' 


Ci, 1 

a n) 


= 


, Ci = 



. fm — 1 _ 


. c n _ 


. Ci ' k i . 


(4.5) 

(* = l,...,n), (4.6) 


where V^(a i, 


i CX- r 


is defined in (11.131) . 
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Lemma 4.2. Given spherical chains (11.161) and associated polynomials (13.61) . the fol¬ 
lowing are equivalent: 

(1) V£(aci,...,at n )X m = 0, where X m = [g 0 ... g m -i] T . 

(2) The polynomial g(z) = gj satisfies homogeneous interpolation condi¬ 

tions 

g\t = 0 for i = 1,..., n; j = 1,..., fc*. (4.7) 

(3) g is in the right ideal (G) r generated by G = lrcm(P 0 , 1 ,... ,P an )- 

Proof: The equivalence (1) (2) is the homogeneous version of the equivalence 

(I4.5l) <t» (|4.6p above. Furthermore, according to (13.31) . conditions (|4.7|) are equivalent to 

<7 = 0 (mod r P ai ) for i = l,...,n, 

where the polynomials P ai are defined in (13.61) . The latter conditions are in turn 
equivalent to the single condition g = 0 (mod r G ). □ 

Corollary 4.3. Given spherical chains (11.161) . let k be the integer defined in The¬ 
orem HZ3 (equivalently, k = deg(lrcm(P ai ,..., P a „)))- Then k leftmost columns in 
V^(aq ,... ,a n ) (m > k) are right linearly independent. 

Proof: Since g = 0 is the only polynomial in the ideal (lrcm(.P Q , 1 ,..., P an ))r of degree 
less than k, it follows by the equivalence (1) (3) in Lemma f4.2l that the homogeneous 

equation (aq,..., a n )X K = 0 has only trivial solution. Therefore the columns in the 
matrix V*(ac i,... , a n ) (i.e., k leftmost columns in VJ~(ai,... , a n ) for m > k) are right 
linearly independent. □ 


Theorem 4.4. Assume that the leftmost elements ocii,...,a n i in given spherical 
chains (11.161) are all distinct and none three of them belong to the same conjugacy 
class. Then the square matrix Vfr(a .\,..., a n ) (N = ]Th ki) is invertible and all solu¬ 
tions to the problem (14.51) are given by the formula 

f = / + Gh, where G := lrcm{P ai ,... , P an }, (4.8) 


where P ai is the polynomial defined in (El]). Ii £ B[z] is a free parameter, and where 


f(z)=[l z ... z N x ] 7jv(ai,...,Q„) 


-l 


Cl 


c n 


is a (unique) solution to the problem. (|4.5I) of degree less than N. 


(4.9) 


Proof: By R emark l3.61 deg(G) = N and hence, the columns of the matrix V^(ol\, ..., a n ) 
are right linearly independent (by Corollary 14.31) . so that the matrix is invertible. We 
now can solve the non-homogeneous equation (14.6|) (with m = N) to get the column of 
coefficients /o, • • •, /jv which being substituted into the formula 

n- i r f 0 - 

f( z ) = J2z j fi = [i ••• ^ _1 ] : , 

j =o L/jv-i. 

leads us to (14.91) . Since Gh is the general solution to the homogeneous problem (14.71) . 
the formula (14.81) follows. Since deg(/) < N (by ()4.9j) ) and deg(G) = N, it follows that 
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deg(/ + Gh) = deg(G) + deg (h) > N for any h ^ 0, and therefore, / is indeed a unique 
solution to the problem (14.51) of degree less than N. □ 


Corollary 4.5. The square confluent Vandermonde matrix V)~ (a i,..., a n ) based on n 
spherical chains (11.161) is invertible if and only if the leftmost elements 01 , 1 , 02 , 1 , • • •, o n ,i 
are all distinct and none three of them belong to the same conjugacy class. 


Proof: The “if’ part has been proven in Theorem 
let us consider the usual Vandermonde matrix V* 


I4~4l 


most elements 01 , 1 , 02 , 1 ,..., a n ,i of the given spherical chains. If the confluent matrix 
VJ~(oq ,... ,a n ) is invertible, the rows in are left linearly independent, so that V ,^ 
is invertible as well. Then the claim follows from Theorem 01 □ 


The formulation of the right Lagrange-Hermite problem and the “right” versions of 
Theorem 14.41 and Corollary 14.51 will be omitted. 


5. Representation formulas for divided differences 

Formulas (11.111) . (11.121) (termed “representation formulas” in [ 6 ]) express the value 
of a polynomial at any point in the conjugacy class in terms of its values at any two 
other points from the same class. It is clear that the same formulas hold for formal 
derivatives of /, that is, for divided differences of / based on the spherical chains of the 
special form (12.81) . Theorem 12.41 asserts that the divided differences of / based on two 
spherical chains of the same length k from the same conjugacy class Sci and with 
distinct leftmost entries define all divided differences of / of order up to k based on any 
chain 0:3 C S. A more general question is: given arbitrary spherical chains aq, ol 2 C S, 
for which spherical chains 0.3 C S does the formula (H3D hold with some matrices A 
and B independent of f ? 


Theorem 5.1. Let cxi = (a*,i,..., a^) (i = 1,2, k 2 < k\) be two spherical chains in 


the same conjugacy class ScB sharing 1 x 2 leftmost elements: 

ai,j — c> 2 ,y (I j T V 2 ) and ex i,^ 2 +i 7 ^ ^ 2 ,u 2 +i- 

Let 0:3 = (0:3,1,..., a 3,^3) be another spherical chain in S such that 

ks < k\ and k^ — 1 x 3 < k 2 — 1 x 2 , (5-2) 

where v 3 is the number of leftmost elements shared by a 1 and 0:3. Then there exist 
matrices A € HI fc3Xfcl and B € H fc3Xfc2 such that 

A^(a 3 ; g) = M f (ai; g) + BA £ (a 2 -, g) (5.3) 

for any polynomial g € El [2]. Consequently, 

V m( a s) = Av m( a i) + BV m( a 2 ) for m = 1, , 2 , .... (5.4) 


Proof: We start with the Lagrange-Hermite interpolation problem 
[ a i,i, • • •, oi,^, f\t r l .j for j 1 ,..., />'j. 

[o 2 ,i, • •. f]i = c 2}j for j = v 2 + l,...,k 2 
based on the chains aq, a 2 ', the target values qj will be specified later. 


(5.5) 
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Step 1 : The polynomial f(z) = z^ fj satisfies conditions (15.51) if and only if 

j =o 


h,i 0 0 

0 0 — V2 


V m (a.i,cx 2 )X m — 


C\ 

C 2 


(5.6) 


where V^cxi^cxfi) is the confluent Vandermonde matrix based on ai, a 2 , 
Ci = 


Cl,l 


c 2,u 2 +l 


' /o ' 


, C 2 = 


and X m = 




. C2 - fc 2 


fm—l 


Indeed, making use of notation (11.131) . we write conditions (15.51) in the vector form as 
A^(ai; /) = Ci, [0 4 2 -i/ 2 ] A^(a 2 ; /) = C 2 , 
or in the matrix form, as 

(5-7) 


4i 

0 0 

Xl(OLl] f) 


c{ 

0 

0 Ik 2 — V2 

_A i(ol 2 \ f) 


C 2 


and therefore, in the form (15.61) . due to (12.181) . 

Step 2 : The square matrix 

V^ 1+k2 _ U2 (a.i, a 2 ) is invertible. 


K = 


I kl 0 0 

o o Ik 2 -V2 


(5.8) 


If we set Cij = 0 everywhere in (15.51) , the obtained homogeneous conditions are equiv¬ 
alent (due to assumptions (J5.1D and the equivalence (13.31) 1 to 


/ = 0(mod r P Q!i ) where P ai = p ail 



1 , 2 ). 


Equivalently, / = 0(mod r G) where G = lrcn^Po^, P a2 )i and it follows from (15.11) by 
Lemma 13.31 that 


G = 


■ P 


— fc 2 +i / 2 


(5.9) 


Since deg(Ts) = 2 , we have deg(G') = ki + k 2 — v 2 and hence / = 0 is the only solution 
to the homogeneous problem (15.51) of degree less than ki + k 2 — v 2 . Equivalently, the 
homogeneous matrix equation 


Iki 

0 


0 

Ik2—V2 


^kl+k2— V2 (**1 ’ OL ‘2)X kl +k2— V2 0 


has only trivial solution. Therefore, the matrix K of this system is invertible. 


Step 3: If f £ H[,s] satisfies homogeneous conditions (15.5p . then it also satisfies 
A^(q: 3 ; /) = 0 for any spherical chain a 3 subject to (15.21) . 

Due to assumptions (15.21) . the polynomial P a3 = p aa 1 ■ ■ ■ p aa fc ^ is a left divisor of 
the polynomial G given in (I5.9p . Indeed, since Xg has real coefficients, it commutes 
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with any polynomial in H^]. On the other hand, Xg = p a pa for each a £ S. Since 
k\ — k 2 + V 2 > k 3 — k -2 + v 2 > z/ 3 , we have 

— yk 2 -v 2 ,, „ 

CT — ^ aijl Pai ;fcl _fc 2+V2 


Pa 1,1 

. . vk 2 -v 2 

* ^1,^3 S ' Ctl, 1^3+1 * ^1 — fc2 +^2 


~Pa 3i i 

Pa 3 ,v 3 Pa 3 ,„ 3 +i ' ' ■ P “ 3 ,/j 3 Pa 3 ,k 3 ' ' ' ^« 3 ,1/3+1 X 



yk 2 -v 2 -k 3 +v 2 

S ra i,„ 3 +i c'ai,fc 1 -fc 2 +„ 2 


_p y^2—V2—k$+V2 n_ . . . n 

S' ' ' a 3 ,fc 3 *^3,1/3 + ! ' a l, 1/3 + I ' OL \,k^_ — fc 2 +i / 2 * 

(5.10) 


If / € H[z] satisfies homogeneous conditions (15.51) . then / belongs to (G) r , by (the proof 
of) Step 2. Then / e {Pa 3 )r (by (I5.10P ) which is equivalent to equality A^(a: 3 ; /) = 0 
(by@). 


We now complete the proof of the theorem. Given a polynomial g € H[z], we let 
l l .j [oq,i, • • •, oq 5 j, foi j 1,..., /ti, 

C2j = [a 2 ,i,...,o:2j; 5k for j = v 2 + 1,... ,k 2 , 

and consider the interpolation problem (15.51) based on this data. By Step 1, the problem 
is equivalent to the matrix equation (15.61) with 


(5.11) 


cv 


C 2 



hi 0 
0 0 


0 

h 2 —u 2 


A*(aq; g) 
A g) 


For m = k\ + k 2 — z/ 2 , this equation takes the form 

hi 0 0 


KX, 


ki+k 2 —u 2 


0 0 Ik 2 - 


2 — 1 / 2 . 


A*(oq; g) 

A^(a 2 ; 9) 


where K is given in (15.81) . Since K is invertible, it follows by Step 1 that the polynomial 


hi 0 0 

0 0 h 2 —u 2 


A*(aq; g) 
Ae(a 2 ; 9 ) 


(5.12) 


f(z) = [1 z ... 3 *i+fe-**-i] K~ l 

satisfies conditions (15.51) . Due to the current choice (15. lip of c t j, we have 

A*(oq; /) = A*(aq; g) and A^(a 2 ; f) = A e (a 2 ; g). 

Hence, the polynomial f = f — g satisfies homogeneous conditions (15.51) and hence, by 
Step 3, 

Ai(a 3 ; /) = A £ (a 3 ; /) - A £ (a 3 ; g) = 0. (5.13) 

We now have from (15.131) . (15.121) . (12.121) and (12.151) . 

Ai{a 3 \g) =At(a 3 -J) 

hi 
0 


[ hf : ... 

=v k ‘ I+la _„ 2 ( a3 )i r- 1 

which implies (15.31) with 

A = H7 +fc2 _, 2 (a 3 )A- 1 


J'^ +fe 2 -" 2 - 1 A fe3 ] A" 


-i 


0 0 
b h 2 —u 2 


A*(oq; 5 ) 
A*(a 2 ; 5 ) 


h 1 0 0 

0 0 Ik 2 —u 2 


A*(aq; g) 
Ae(a 2 ] g) 


hi 

0 


B = v£ l+k 2 _Jcx 3 )K 


-1 


0 0 

b h 2 —i/ 2 








































19 


Since formula 115.31) holds for any polynomial g G H[*], we have, in particular, 

A^(q: 3 ; z j ) = z j ) + BA e (a 2 -, z j ) for all j > 0 (5-14) 

and, since A^oq; z j ), Ai(a 2 -, z J ), A^(aq; z J ) are the j-th columns in the matrices 
V^(aq), Vm(a 2 ), V^a^) respectively, equality (15.41) follows. □ 

In terms of Section 3, the last theorem can be reformulated as follows. 


Remark 5.2. Given two spherical chains oq, a 2 in the same conjugacy class Sci, let 
P a i and P a2 be the associated indecomposable polynomials (|3.6D . Then for any spher¬ 
ical chain a 3 whose associated polynomial P a3 is a left divisor of the lrcm(P Q , 1 , P a2 ), 
equality (15.31) holds for any / G and matrices A, B independent of /. 

Theorem 11.41 is a particular case of Theorem 15.11 for which, however, we will present 
more explicit formulas for the coefficient matrices A and B. With a spherical chain a, 
we associate the square matrices 

V a := V k f (a) and T a = V~ l J k V a , a = (oq,. .., a k ), (5.15) 

where J a is given in (12.121) . Recall that V a is invertible as the square upper triangular 
matrix with all diagonal entries equal one. 


Theorem 5.3. Under the assumptions of Theorem \l.f\ equality (11.141) holds for every 
f £W[z] with 

,- V -Y (5-16) 


A — Va 3 (Ta. 3 T a2 )(T ai T a2 j , ai , 


B = V a3 (T a3 - T ai )(T a2 - T a 


OL2 ’ 


where V ai and T ai are k x k matrices defined as in (15.151) . 


Proof: We get Theorem 11.41 bv letting k\ = k 2 = k 3 = k and v 2 = 0 in Theorem 15.11 
Therefore, the formula (11.141) holds with 


A = V&(a 3 )K 


-i 


B = Vf k (a 3 )K 


-l 


(5.17) 


where, according to (15.81) . K = V^,(ai 1 , 0 : 2 )- Making use of notation (12.151) . we can 
write 

4(« 3 )=Y„3 J^V a3 \=V a3 [I k T a3 ] (5.18) 

and similarly, 

V£ k {a.i,a 2 ) = 


1 

0 

A* 


4 T ai 

0 Vck 2 


4 T a2 


Since a 2 ) is invertible, the right factor on the right side of the latter equality 

is invertible and hence, (T a2 — T ai ) is invertible. We have 

\T a2 -T ai )-'V ai l 0 

0 (T Q2 - T a 


Ta. 2 

-T ai ' 


-4 

Ik . 



Vf k (a u a 2 ) 1 = 

which together with () 5.18 1) leads us to 
v 2k ( a 3) V2k ( Q i > Q^) -1 


-OLi) Va . 2 J 


T a 2 

1 

_ 1 


-4 

4 



(T a2 -T ai )-'V-? 

0 


y “3 \- Ik r « 3 J [-4 4 J [ 0 (T Q2 - T ( 

= v a3 [(T aa - T a3 )(T a2 - T a f)-'V^ (T a3 - T ai )(T a2 - T ai )~ l V~ 2 ] . 


Oi-1 ) V OL2 J 
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Substituting the latter equality into (I5.17p . implies (15.161) . thus completing the proof 
of the theorem. □ 

Example 5.4. Let k = 1 and Qj = (a*) C S for i = 1,2,3. According to (15.151) 
and (12.121) . V ai = 1 and T ai = J ai = a*, while formulas (11.13|) and (12.51) show that 
A i(cxu f) = [oti', f]e = f ei (ai). In this case, formula (j!.14jl amounts to (11.111) . 

It can be shown that the matrices (I5.16P are lower triangular (although the matrices 
V ai is upper triangular and T ai is not triangular at all). Furthermore, the spherical 
chains an = (a,..., a) and a .2 = (a ,... ,a) of the form (11.31) satisfy the assumptions 
of Theorem 11.41 and, on account of (12.81) . we arrive at the following conclusion. 

Remark 5.5. The left divided difference [ 71 ,... , 7 y, f]i based on the spherical chain 
7 = ( 7 i,..., 7 fc) C S' is a left linear combination of 2k elements (f^) ee (a) and 
(/ W ) e * (a) (j = 0 ,..., k — 1 ) with the coefficients depending on 7 and a, where a 
is an arbitrary fixed element in S. 


Remark 5.6. If 71 = ... = 7i , = a ± 7 ^+ 1 , then [ 71 ,... 77 f]t equals for 

j = 1 ,..., v (by (12.81) 1 and it is a left linear combination of (f^) ee (a) (j = 0 ,..., j — 1 ) 
and (/^) e ^(a) (i = 0 ,..., j — v — 1 ) for j > v. 


Example 5.7. Let 07 = ( 07 , 07 ), a .2 = ( 02 , 02 ) and 0:3 = ( 01 , 02 ) for two elements 
«i ~ 02 £ H. Then 


Foil - hc«3 — 


T = I / -1 T 2 V = 

v n/o y olz 


l «3 


0:3 ^ OL3 v a 3 


and similarly, T ai = 


1 Ol 
0 1 

1 —Ol 

0 1 

—a 2 —2of 

2oi 3o? 


V aa = 


1 02 

0 1 



p 

1 —‘to 

, 0 


1 Ol 


— Q7CI2 —Oio| — 0l0 2 0l 


OL 2 + Ol OL2 


0 1 


02 + «1 o| + a 2«l + 


T — 

± OL 2 - 


Oia 

2 q 2 


2 - 2 a?' 


3«9 


By (12.81) . the equality ()1.141) now takes the form 


/ e *(« 1 ) ' 
[ 01 , 02 ; f]t 


= A 


+ B 


/^(o 2 ) 

(/O^M 


(5.19) 


' f ei (a 1 ) 

_(m« 1 ) 

where (as quite tedious calculations of the right hand side expressions in (15.161) show) 

0 


A = 


B = 


(o 2 — 0 2 ) 


-1 


1 0 
(o 2 o 2 ) —1 (o 2 - a 2 )^ 1 (o 2 - Ol) -1 

It is readily seen that the equality of the top entries in (15.191) is trivial, while the 
comparison the bottom entries proves formula ( 12 . 101 ) . 


Proof of Theorem 11.71 Given spherical chains (11.161) . let k be the integer defined in 
Theorem 11.71 By Corollary 14.31 the n leftmost columns in the matrix V* l (a \, • • •, a n ) 
(for m > k) are right linearly independent. Thus, rankV^(a 7 ,... ,a n ) = m if m < n 
and rank (07,... , a n ) > k if m > k. It remains to show that 

rankV^(ai, ..., a n ) < k if m > k. 


(5.20) 
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Let S C HI be a conjugacy class that contains at least three spherical chains from (11.161) . 
Let Q.i be a chain with the maximal length and for any other chain ctj C S, denote by 
i/j the number of leftmost elements shared by ctj and aq. Let ot r = (a rj i,..., ot r ,k T ) be 
the chain in S for which the integer kj — Vj is maximally possible. In terms of associated 
polynomials (13.61) . we choose two polynomials P ai and P ar with maximally possible 
degree of their lrcm. 

Any spherical chain ctj C S different from cq and a r satisfies conditions 

kj < ki and kj — Vj < k r — v r . 

By Theorem 15.11 equality (|5.4I) holds for all m > 1. So, removing the block Vf n (ctj) 
from V^(aq,... ,ot n ) does not change the rank of the latter matrix. Repeating this 
argument, we remove from V^(ai,..., a n ) all blocks V^JjOLj) based on the chains in S 
different from ctj and a r . We then observe that the v r top rows in V^joti) are identical 
to the corresponding rows in V^(a r ). Hence, removing the u r top rows from the block 
V^(o: r ) does not change the rank of the matrix. After this removal, the remaining 
matrix is of the same rank as the original one but contains only two blocks (V^(aq) 
and [0 Ik r ~v r \ ^m( Q r)) associated with the elements in the class S. The total numbers 
of rows in these blocks equal (by Proposition 13.51) 

ki + k r — v r = n(S) = deg(lrcm(P aj : ctj C S)). 

Repeating the latter procedure for each conjugacy class Sj containing more than two 
spherical chains from (I1.16p . we come up with the matrix of the the same rank as the 
original matrix V^(ai,... ,ot n ) but having only k = ^ ■ n(Sj) rows. Its rank cannot 
exceed k which implies (12.211) and completes the proof of Theorem 11.71 □ 


6. Formal power series over quaternions 

In this section, we discuss several analogs of the preceding results in the context of 
the space H[[z]] of formal power series over H. Given an / E H[[~]], we denote by f* is 
conjugate power series defined by 

oo oo 

f\z) = Y J ^j if f( z ) = Yl z3 fr l 6 - 1 ) 

j =o j=0 

The anti-linear involution / i—>■ f* can be viewed as an extension of the quaternionic 
conjugation aoa from HI to HI[[z]]. 

6.1. Linear independence of certain power series. We consider the power series 

OO OO 

k a (z) = ^ ~^a k z k and k^(^) = k.a(z) = ^ ^a k z k (a € HI). (6.2) 

k=0 k =0 

Given a spherical chain a = (aq,..., «&), we define the infinite matrix V^ 0 (a) by letting 
m = oo in (1 1.5 [) and we use the entries from the same row in V£,(ot) to define power 
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series 


oo oo 

fi( z ) = Y^ [«1 '■’ zJ ]e zJ = = 

j =0 j=0 


f 2 (z) = y [a 1 ,a 2 \z 3 ] t z 3 

3=0 



zk a2 (z)\t ai (z), 


(6.3) 


fk(z) =Y, [ai, • • .,a k -,z 3 ] t z 3 = z k 1 k afc ( 2 :) • • • k a2 (z)k Ql (z). 
j=o 

By (16.11) . (16.21) and (12.191) . the conjugate power series /■ are given by 

fj{z) = z 3 - 1 ^ (z)--- kg,. (z) = z?- 1 ^ (z)--- k s . (z) (6.4) 

OO OO 

= y [«1. ■ ■ -^j-z^z 1 = y [ail, ■ ■ ■ ,aj-,z l ] r z l (j = 1 
1=0 2—0 

Proposition 6.1. Given spherical chains cti = (a^i,..., a * >ki ) (i = 1,..., n), let k be 
the integer defined in Theorem m Then the dimension of the left linear span of power 
series 

fi,j(z ) = z 3 ~ 1 k a . J (z) • • k ai l (z) (i = 1,... ,n; j = 1,..., k { ) (6.5) 

and the dimension of the right linear span of the conjugate power series 

/y'( 2: ) = ^' lk 5 i,iW'" k 5 il3 '(z) (i = l,...,n; j = l,...,ki) (6.6) 

equal k. In particular, the series (16.51) are left linearly independent (the series (|6.6[) 
are right linearly independent) if and only if the leading elements OL\ y i,... ,a n> i are all 
distinct and none three of them belong to the same conjugacy class. 

Proof: By (16.31) . the coefficients of the series (16.51) are the entries from the same 
row in the infinite confluent Vandermonde matrix V^(a i,..., cx n ), and the statements 
concerning the series (16.51) follows from Theorem H771 The statements concerning the 
conjugate power series are now immediate. □ 

6.2. Quaternion formal power series and their evaluations. Most of the pre¬ 
ceding results invoked left and right evaluations of quaternion polynomials. To extend 
them to the setting of H[[z]], we restrict our attention to a class of power series for which 
left and right evaluation functionals make sense. We denote by Mr = {a £ H : |a| < R} 
the open ball in H of radius R centered at the origin, and we introduce the space 

Hr = j/(v> = y f] z3 ■ linrsup f/\h\ < j. 

Remark 6.2. The power series k Q (16.21) belongs TLr with R = |a| _1 . More generally, 
the power series z 3 1 k Qj . (z) ■ ■ • k ai (z) belongs to Hr with R = min{\ai\ 1 }. 
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We also note that for any f(z ) = fkZ k in Hr and any a £ Br, the series YlT=o ak fk 
and Y^k =o fkOi k converge absolutely, so the evaluation formulas (11.511 (with m = oo) 
make sense. Furthermore, the power series L a f and R a f (defined as in (11.61) but with 
m = oo) are also in Hr. Therefore, left and right divided differences for / £ Hr can 
be defined via formulas (i23D. m for any elements a\,... ,a k £ Br. 


Remark 6.3. Representation formulas for quaternion polynomials in Theorems \l.f\ 
fO and hold true for all elements g £ Hr and spherical chains q 1 ,q: 2) Q: 3 €= Mr. 

Proof: We verify (the most general) Theorem 15. 11 For g(z) = gjZ J we have from 

(15.141) by linearity, 


A^(a 3 ; 9) = Y A ^ a3 ' z ^9j = ^(^A £ (ai; z J ) + BA e (a 2 ; z 3 ))gj 
1=o 1=o 

= AAi(ai; g) + BAi(a. 2 ] g ), 

where convergence of all series follows since au, Q 2 , 0:3 £ Br. 


□ 


6.3. Square summable formal power series. We now consider the space 


H 2 = \ h(z) = zjh j '■ 
j =0 


\h 


Y 1 ^' 


< 00 


J=0 


of elements in H[[z]] with square summable coefficients, endowed with the left and right 
inner products 

OO OO 

(h, g)t = Y^j h T (■ h . 9 )r = Y h ^r 

j =0 i=o 

It is clear that Hr C H 2 C H\ (for each R > 1) as sets. The series (16.21) belongs to H 2 
if and only if a £ Bj (i.e., \a\ < 1), and is of particular interest due to the following 
reproducing property: 

OO 

(h, = Y* k h k = h ei (a) for all h £ H 2 . 

l=o 

More generally, if a\,,aj belong to Bi, then the power series 

fj(z ) = z 3 ~ l k aj (z) ■ ■ ■ k ai (z) 

belongs H 2 (by Remark 16.31) and reproduces the j-th divided difference 


(h, fj)i = otj‘,h] t for all h £ H 2 , 

which is verified by a straightforward power series computation. 


( 6 . 8 ) 


6.4. Cauchy matrices. Given spherical chains cv, = (ccj i,..., a^) C Bi (i = 
1,..., n), the associated the power series (J6.5I) and their conjugates (| 6 . 6 D belong to H 2 . 
Let us denote by Pj- the left Gram matrix of the set 


T = {fL( z ) : * = L, • • •, n; j = 1, • • •, ki} 


( 6 . 9 ) 
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where ffjiz) = z ] 1 k Si>1 (^) • • • k a u {z). Thus, 


Pr = 


= Where 


P ~ 

1 OLi.OL-f — 


f t \ y- 1 ’~s 


( 6 . 10 ) 


Remark 6.4. In case ki = 1 for i = 1 we have T = {k Sl ,..., k Sn }, and 

consequently, 


P =[(k 5 ,.,k a , )*]”,,= £ 


<x\a) 


L k=0 


( 6 . 11 ) 


i,j =1 


If ai,aj = ctjcti, then Ylh=o a i®j = (1 ~ adT/) -1 - For this reason, we will refer to the 
matrix (16.111) as to left Cauchy matrix and we will call the matrix (|6.10l) a generalized 
(or confluent) left Cauchy matrix. 

Connections between quaternion Cauchy and Vandermonde matrices were observed 
in [1]. The next theorem summarizes the confluent case. 


Theorem 6.5. Let Pj be the confluent Cauchy matrix based on n spherical chains 
on = (a^i,..., otiki) C Bi (i = 1 and defined as in (I6.10p . Then Pj, can be 

factored as 

P T = V£(ai, • ■ ■, a n )V^(a i,..., a„)*, (6.12) 

where V^,(ai,... ,ot n ) is the infinite confluent Vandermonde matrix based on the same 
chains. Therefore, 

(1) rankPjr = rankV^(aq,..., a n ). 

(2) Pjr is positive definite if and only if the leftmost elements aqi,... ,a n ,i are all 
distinct and none three of them belong to the same conjugacy class. 

(3) The block Pi p of Pj satisfies the Stein equality 

P oti,oi l i ~ T Q , i P Q , iiQ , i ,X a . / = Ep.E k / , (6.13) 

where Z ai , Z a .,, E^., E^., are defined via formulas (|2.12p . 


Proof: Due to the block structure (16.101) and (11.171) of Pjr and V^(qi, ..., a n ), in 


■ ini ui P i 

order to prove (| 6 . 12 p it suffices to verify equalities 

Pai,a it = V<L( a i) v 2( a i'Y for i,j = 1 ,... ,n. 
To this end we first observe that on account of formula (1 1.15 f) . 


(6.14) 


Fqo (^i)koo ^ ^ [th .1 1 ■ ■ ■ i &i,j i z ^ i , . . . , Oip j' , z ] «. (6.15) 

L s =0 

On the other hand, by virtue of formula (15.21) . 

oo oo 

fife) = Y = Y [ ai ’p ■ ■ ■ z8 ]i z8 ’ 

s=0 s=0 

which imply, by the definition (16.71) of the left inner product in H 2 , 

oo 

= I = Yh ^ ai di ■ ■ ■ > a ij'i zS \l ' \ a i', 1) ■ • • ) a i\j r i zS ](_■ 


s =0 


Comparing the latter equality with ()6.15l) completes the verification of (16.141) . 
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Statements (1) and (2) follow from factorization (16.121) and Theorem 11.71 To prove 
()6.13j) . we start with the Stein identity (12.171) (with m = oo) applied to the chains on 
and OLi> : 


V&tXi) = J^V^Fl + E ki E^, *£(<**) = J a Mcx t ,)F^ + E k ,E 


rT 


i 

oo v 


ri 

' oo v 


oo 


pT 


Multiplying the first equality by the adjoint of the second equality on the right gives 


= (j ai Vi(ai)FZ + E ki El) (F^V^cy,)* , 
which, on account of equalities 


f t f t — t f t f — n f t f — ^ 

r oo r oo ~ 1 J r oo-^oo ~ u, -E' 0 o- C/ oo — 


simplifies to 


V^aJV^cXi')* = +E kt Ej, 


The latter equality is equivalent to (J6.13D . due to factorization (I6.14p . □. 

The right generalized Cauchy matrix Pf, is defined by the formula (16.101) but with the 
right inner product rather the left. The equality Py = V^(ai,..., ctn^V^au,..., a n ) 
is verified along the lines of the proof of Theorem 16.51 


6.5. Lagrange-Hermite interpolation. To impose the Lagrange-Hermite interpola¬ 
tion conditions (14.5[i we need well defined evaluation functionals. If the given spherical 
chains (11.161) are all in Mr, we may formulate the interpolation problem (14.51) for power 
series from the space 'Hr. Since the formulas (11.41) (and therefore, more general for¬ 
mulas (12.41) 1 hold for any / € Hr and a\,...,a k € Mr (see e.g., [BJ[8]), it follows that 
any power series g £ Hr satisfying homogeneous conditions (14.7|1 is necessarily of the 
form g = G ■ q for some q £ Hr where G is the polynomial defined in (14.81) . Since the 
problem (14.5|1 is linear, we arrive at the following result. 

Proposition 6.6. Given spherical chains (11.161) in B/?, such that their leftmost ele¬ 
ments aiq,..., a n ,i are all distinct and none three of them belong to the same conjugacy 
class, a power series f belongs to Hr and satisfies conditions (1431) if and only if 

f = f + Gh for some h&HR, (6.16) 

where f and G are the polynomials defined in (14.91) , (14.81) . If R = 1, then f of the form 
(|6.16l) belongs to H 2 if and only if h £ H 2 (since f and G are both polynomials). 


6.5. Norm-constrained interpolation in H 2 . We finally consider the problem 
of finding all power series / £ H 2 satisfying interpolation conditions (|4.5I) and the 
additional norm constraint ||/||h 2 1- To relate (fairly explicitly) the norms of / and 

the corresponding parameter h, the formula (|6.16l) should be modified. The first step 
is to replace the particular solution / by another one which is orthogonal to G ■ H 2 , 
the solution set of the homogeneous problem (14.71) . To this end, observe that by the 
reproducing property (16.81) . interpolation conditions (|4.5j) can be written as 

(/> f[j)t = (H,j for i = j = 


(6.17) 
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where /' ? are defined in (16.51) . Thus, g G H 2 satisfies homogeneous conditions (14.71) 

(i.e., g G G ■ H 2 ) if and only if it is left-orthogonal to all elements /? ■ from the set (16.91) 
and hence, to their right linear span 

M = span r | ffj ( z ) ■ i = l,...,n; j = 

In other words, Ai 1 = G • H 2 , and thus we are looking for a power series f m \„ G Ai 
satisfying conditions (16.171) . Writing / m ; n in the form 

n ki 

fmin(z) = ^2 fiA Z ) d V ( 6 - 18 ) 

i =1 J=1 

with unknown coefficients dij € H, we take its left inner product against all elements 
from (16.91) getting (on account of interpolation conditions (16.171) 1. the linear system 




dij — 


Ci,' 


1 ,3 


for all % = 1. 


,n]j — 1,..., kv 


(6.19) 


*=i i=l 


The matrix of this system is Pjr which is invertible under assumptions of Proposition 
16.61 Substituting the coefficients di.j found from (16.191) into ()6.181) we get the formula 


fmm(z) = [. Fl(z ) 


F n {z)\ (Pjr) l C, where C 


'CW 

.C n _ 


( 6 . 20 ) 


is the column defined in (14.61) and where 


Hz) 



for i = 1,..., n. 


It is not hard to show that for f m - m of the form ()6.20j) . || n || 2 p = C*(Pjr ) l C. 

The next modification concerns the second term on the right side of (|6.16j) . that is, 
the general H 2 -solution of the homogeneous problem (14.7p . 


Proposition 6.7. H 2 = k/j • H 2 (as sets) for any fi € and kg defined as in (16.2p . 

Proof: For any h G H 2 , the series g(z) = (z — fi)h(z) also belongs to h G H 2 and, since 
k q(z)(z — (3) = 1, we have h = k pg so that h G k^ • H 2 . Thus, H 2 C k^ • H 2 . 

To verify the reverse inclusion, take h(z) = Y^k zk G H 2 and let h(z) = Y \^k\z k ■ 
Then h G !R[[2r]] actually belongs to the Hardy space H 2 (0) of the unit disk of C and 
obviously satisfies = ||/i || H 2 = II^-IIh 2 < oo. Since the power series \a.\p\(z) = 

Y \fi\ k z k is a bounded in the closed unit disk, 

2 oo , k x 2 

— XT ( ‘ \dj \ J 

k =0 'j =0 2 

= ll k |/3|^ll|/ 2 (ro) < (™l k |/3|l) 2 • IMHfio) = (1 - I/5|)- 2 ||^IIh- 

Therefore, k gh belongs to H 2 verifying the desired inclusion and completing the proof. 

□ 


IMS, = £ 




















27 


Proposition 6.8. For any monic polynomial G E H[^] of degree d > 0 and having no 
zeros outside Bi, there exist the elements (3\, ..., (3d (each (3j is equivalent to one of the 
zeros of G) such that the power series 

0 = G ■ k i3l ■ k 02 ■ ■ ■ k 3d (6.21) 

has the following property: || 0 / i || H 2 = ||/i || H 2 for all h E H 2 . 

We refer to [U Theorem 7.1] for the proof and explicit construction of (3\,..., (3d- In 
the commutative case, 0 is just the Blaschke product having the same zeros (counted 
with multiplicities) as G. It follows from Proposition 16.71 and formula (16.211) that 
G ■ H 2 = 0 • H 2 and hence, G can be replaced by 0 in the parametrization formula 
(16.161) . We thus arrive at the modified parametrization formula 

/ = /min + 0!i for some h E H 2 (6.22) 

with fmin given by (16.201) and 0 constructed as above. The advantages of the formula 
(16.221) can be seen from our last theorem. 

Theorem 6.9. Under the assumption of Provosition [6i~6l the formula (|6.22l) describes 
all f E H 2 satisfying interpolation conditions (14.51) . Moreover, the representation (16.22[) 
is left-orthogonal and therefore, 

II/IIh 2 = 11 /min 11H 2 + II^IIh 2 = ^(P^C + \\h\\ 2 ^. (6.23) 

In particular, fmin is a (unique solution) to the problem (14.51) with the minimally possible 
norm. In case C*(P^)~ 1 C < 1, all solutions to the problem with ||/|| H 2 < 1 are given 
by formula (|6.22l) with h subject to ||/i||| 2 < 1 — C*(Pj r )~ 1 C. 
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